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ABSTRACT 



Context. According to the modem cosmological paradigm, cosmic voids form in low density regions between filaments of galaxies 
and superclusters. 

Aims. Our goal is to see how density waves of different scale combine to form voids between galaxy systems of various scales. 
Methods. We perform numerical simulations of structure formation in cubes of size 100, and 256 h~ l Mpc, with resolutions 256 3 and 
512 3 particles and cells. To understand the role of density perturbations of various scale, we cut power spectra on scales from 8 to 
128 /T 1 Mpc, using otherwise in all cases identical initial random realisations. 

Results. We find that small haloes and short filaments form all over the simulation box, if perturbations only on scales as large as 
8 lr l Mpc are present. We define density waves of scale > 64 h~ l Mpc as large, waves of scale = 32 h~ l Mpc as medium scale, and 
waves of scale ^ 8 h~ l Mpc as small scale, within a factor of two. Voids form in regions where medium- and large-scale density 
perturbations combine in negative parts of the waves because of the synchronisation of phases of medium- and large-scale density 
perturbations. In voids, the growth of potential haloes (formed in the absence of large-scale perturbations) is suppressed by the com- 
bined negative sections of medium- and large-scale density perturbations, so that their densities are less than the mean density, and 
thus during the evolution their densities do not increase. 

Conclusions. The phenomenon of large multi-scale voids in the cosmic web requires the presence of an extended spectrum of primor- 
dial density perturbations. The void phenomenon is due to the action of two processes: the synchronisation of density perturbations of 
medium and large scales, and the suppression of galaxy formation in low-density regions by the combined action of negative sections 
of medium- and large-scale density perturbations. 

Key words, large-scale structure of the Universe; early Universe; cosmology: theory; methods: numerical 



1. Introduction 

The goal of this series of papers is to study the role of per- 
turbations on vario us scales in formation the cosmic web. 
lEinasto et alj ([20 11) used wavelet techniques to understand the 
formation of rich systems of galaxies - clusters and superclus- 
ters. They conclude that superclusters are objects where density 
waves of medium and large scales combine in similar phases to 
generate high density peaks. Similarly, voids are regions in space 
where medium- and large-scale density perturbations combine in 
similar under- density phases of waves. Su hhonenko et alJd201 lb 
demonstrated that the properties of the cosmic web depend es- 
sentially on density perturbations of small and medium scales, 
whereas perturbations of large scale > 100 h^ 1 Mpc modulate 
the richness of galaxy systems from clusters to superclusters, 
and make voids emptier. This paper is devoted to the study of 
the influence of medium- and large-scale density waves on the 
structure and evolution of voids in the cosmic web. 

The cosmic web was first openly discussed at the 
IAU Symposium on Lar ge Scale Structure of the Universe 
dLong air & Ei nastoL Il978l) . At this symposium, four groups re- 
ported results of studies of the three-dimensional distribution 
of galaxies in space using available data for the redshifts of 
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galaxies. The presen ce of voids in the distribution of galax- 
ies wa s reported by Joeveer & Einasto! (119781). iTarenghi et alj 
(Il978h . Ffifft & Gregory! (1 19781) . and iTullv & Fisher! (1 19781) in 
the Perseus-Pisces, Hercules, Coma, and Local superclusters, re- 
spectively. 

The main results reported in this symposium were that: 
(1) galaxies, groups, and clusters of galaxies are not randomly 
distributed but form chains, converging in superclusters; (2) 
the space between galaxy chains contains almost no galaxies 
and forms voids of diameters 20 . . .70 h Mpc; (3) superclus- 
ters are not isolated systems, but are connected by galaxy fil- 
ament s to a connected network - the supercluster-void net- 
work jjjdeveer et all [197! lEinasto et all [T980t IZeldovich et all 
1982; IOortlll983h or the cosmic web (iBond et all 1 1996b . These 
early results were confirmed by the Second Harvard Sky Survey 
( de Lapparentetail 119861; iGeller & Huchrall 1989b and the dis- 
covery of the very large Bootes void by iKirshner etaD(fT98Tl 
Il987h . 

The presence of long essentially one-dimensional galaxy and 
group/cluster chains, and the elongated shape along the chain 
of the central cluster galaxies (often supergiant galaxies of type 
cD) suggests that galaxies and groups/clusters of the chain had 
formed within the chain simultaneously with the formation of 
the whole cosmic web. This occurred in the gaseous phase of 
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the structure evolution, which allowed the dissipation and can- 
celling of velocities perpe ndicular to the chain ax is. These data 
gave strong support to the IZeldovichl (1 19701 1 19781) pancake sce- 
nario of galaxy formation. The observed distribution of galaxies 
was quite similar to the distribution of simulation particles in 
a two-dimensional numerical simulation of the evolution of the 
structure of the Universe, prepared by Shandarin (1975, privat e 
communication) and published by iDoroshkevich et all 0980). 
In this simulation a network of high- and low-density regions 
was seen: high-density regions form cells that surround large 
under-dense regions. S ubsequent three-dimensional simulations 
confir med this picture (iKlvpin & ShandarinlJl983HMelott et al.L 
119831) . 

However, some important differences between the model and 
observations were evident. First of all, there exists a rarefied pop- 
ulation of simulation particles in voids that is absent in real data. 
This was the first indication of physical biasing in galaxy for- 
mation (the term "biased galax y format ion" was introduced a 
few years later bv lKaiserl (11984 . see also lBardeen et all ([1986)). 
A theoretical explanation of the absence of galaxies in voids 
was given by Einas to et al.l dl9 80) (see also a more detailed dis- 
cussion by Einasto et al.l (Il994al) l A simple analytical solution 
of the cosmological evolution of the density of matter shows 
that in over-dense regions the density increases unti l the mat- 
ter col lapses to form compact objects (pancaking by Zeld ovichl 
( 1970)). In contrast, in under-dense regions the density decreases 
but never reaches a zero value - gravity cannot evacuate voids 
completely. 

Even early studies proposed that the cosmic web has a hi- 
erarchical structure. The parts of the web formed by objects of 
different mass or luminosity have different characteristic sizes. 
As voids are defined by the web, the void assembly has several 
important properties: voids defined by more luminous (or mas- 
sive) objects have larger diameters, and voids defined by clus- 
ters/galaxies of certain luminosity contain substructure formed 
by less luminous objects. 

According to the present cosmological paradigm, all 
structural elements of the Universe were formed by the 
growth of initial small density perturbations created dur- 
ing the very early phase of the evolution of the Universe. 
To these elements belong galaxies, groups, clusters, and 
their systems, such as filaments and superclusters. Rich 
superclusters form a cellular distribution, with large voids 
surrounded by rich superclusters. The characteristic diam- 
eter of these superv oi ds is of the or der o f 100 h~ l Mpc 
1 19781 iKirshner et all 1 198 it lEinasto et all 
Su pervoids are not empty, but contain a hi- 
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Superclusters are connected by filaments of galaxies, 
i.e voids have substructure. Observ ationally, this was al- 
ready evident in early v oid studies dJoeveer & Einastol 119781; 
Ide Lapparent et all Il986l) . Theoretica l discussion of t he voi d 
substructure has bee n prese n ted by iRegos & Gellerl dl 989 
Mattel & Wasserman ( 1990), van de Weygaert & van Kampe 
(ll993l) .lGoldberg & Vog elevl(l2004l) :lGoldberg et al. (2005Mnd 
many other s . The s keleton of the cosmic web w a s discussed by 
Hahn et all d2007l). iForero-Romero et al.1 (120091). Sousbi e et all 
(l2008Ll2009l).lArag6n-Calvo et al.l d2010l).lBond et all (l2010allbh . 
IShandarinl d2010l) , and lEinasto et al.l d201 lb . 

It is generally accepted that the initial density perturbations 
had a smooth, extended power spectrum (quasi-flat in terms of 
metric perturbations n s w 1) and a random (Gaussian) distri- 
bution of perturbation phases. The amplitude of perturbations 
(A oc k^Pik)) is larger at short wavelengths and per 6 In k, where 
k is the wavenumber, and P(k) is the power spectrum of pertur- 
bations. For this reason, small objects (mini-haloes in numerical 
simulations and dwarf galaxies in the real Universe) should form 
first. The early formation of dwar f galaxies is con firmed by ob- 
servation of very distant galaxies (Bec kwith et all 120061) . These 
early galaxies grow by the attraction of mor e primordial m atter 
and by clustering, as originally suggested by Pee blesl(ll971l) . 

Initial small-scale perturbations were present everywhere, 
and this raises the question: why do voids not contain galax- 
ies, eve n dwarf galaxie s? This question was asked more specifi- 
cally by Peeblesl d2001l) . The deepest voids in the simulated dark 
matter distribution are never completely empty; they still con- 
tain low-mass condensations of primordial matter - mini-haloes. 
Therefore, voids could be the environment in which faint dwarf 
galaxies are most likely to reside. Why is this not the case? 

This "Peebles question" has stimulated a number of studies, 
both observational and theoretical, to find dwarf void galaxies 
and to study either analytically or by numerical simulations void 
regions of the Univer s e. Studies of void galaxi e s were made by 
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Roias et alJ 00041). iHoeft et al 



Gottlober et al 



Stanonik et al 



iKreckel et alJ d201 lhlah . an d others. Among recen t 
large-scale s u rveys, we mention he r e iHoyle & Voge[eyl (12004 ), 
ICroton et aD d2004l). IConrov et al.1 d2005h. iPatiri et all (I2 006). 
and lTinker et all d2007l) . lKarachentsev et alJ d2003l 120041 12007 ) 
studied the Local Volume out to a distance 10 h' 1 Mpc, 
where they found about 550 mainly very faint galaxies. A 
void galaxy survey (VGS) was initiated by van de Weygaert 
and collaborato rs in our local neighbourhood ou t to re d- 
7 = 0.02 dStanonik et all 12009b Ikreckel et all 1201 lalfht 



shift 



Among the early studies of the void evolution we men 



Ivan de Weygaert et all 1201*1 ). These studies suggest that ex 
tremely faint galaxies do not fill voids, but are located in the 
vicinity of brighter galaxies. 

To investigate the properties of galaxies in voids, the halo 
occupation distribution model is used to determine th e rela- 
tionsh ip between galaxie s and da rk matter haloes (se e JSeliak 
(I2000I). ICoorav & Shethl d2002l). iBerlind & Weinbergl (|2002|) 
Zehavi et all d2004l) . IZehavi et al.1 d2005l) IZheng et al.l (12007 ). 



tion Hoffman & Sha haml dl982l):lHoffman et alJ dl983l).lPeebles| van den Bosch et all d2007l).lTinker et al.1 d2006l 
( 1982). llckddl984 . lFillmore & Goldreichl(ll984l) . lBertschingeri Iwhite et al.1 d2007l) . and IPadmana bhan et al] ( 



d!987l) . among others. The study of the hi e rarchi c al evolution 
of voi ds was pioneered by Dubinski et al.1 dl99 3). Sahni et aD 
( 1994) described the ev olving void hierarchy within the context 
of the adhesion theory. Shet h & van de Wevgaertl (120041) devel- 
oped the excursion set (extended Press-Schechter) description 
of a void hierarchy in the dark matter distribution, followed by 
Furl anetto & Pirar] (120061) comparing this with the void hierar- 
chy in galaxy populations. 



, 2007, 2008a,b), 
and IPadmanabhan et alJ d2009l) ). The re- 
sults obtained by analysing of observational data are in good 
agreement with the results of semi-analytic models, hydro- 



dynamic cosmologic al simulations, and h i gh-resolution colli 

H4|Zheng etall 



signless simu l ations dKravtsov et al. 
IConrov et al 



2004 



2006; von Benda-Beckmann & Mtiller, 



2005 



2008 



Tikhonov & Klvpin, J2009J). In particular, the simu lations of 
Tikhonov & Klypin d2009t) and[lshiyama et"aH ( 1201 ll) show the 
presence of very low-mass haloes in voids, but their mass is 
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probably too low for star formation to be possible. The compari- 
son of the distribution of model haloes and dwarf galaxies shows 
that haloes should have a circular velocity of a t least ~ 35 km/s 
in order to form a galaxy. Tinker et al. (2008b) showed that the 
smallest void haloes have so little mass that no galaxy formation 
is expected according to the model. 

The absence of dwarf galaxies in voids can also be explained 
by gasodynamical p rocesses. Using h igh-r esolution hydrody- 
namical simulations. iHoeft et alj (12006). and lHoeft"& Gottlober 
(1201 Oh showed that photoionisation by the UV radiation field 
is able to stop the cooling and collapse of gas in dwarf galaxy 
haloes. At the redshift z — 0, the characteristic mass scale of 
photo-evaporation corresponds to a circular velocity ~ 27 km/s, 
in good agreement with other studies cited above. 

All studies suggest that galaxy formation is a threshold phe- 
nomenon. In many of the studies cited above, it was demon- 
strated that galaxies do not form in voids, when the void haloes 
have very low masses. However, these studies did not explain 
what is the physical reason for the difference between the typical 
halo masses in the void and supercluster environments? In this 
paper, we try to find an explanation for the difference between 
the masses of haloes in various global environments. 

In the case of the presence of an extended primordial per- 
turbation spectrum, systems of galaxies are produced by an in- 
terplay of density perturbations on all scales. It is clear that the 
difference between the present distribution of galaxies and the 
very early distribution of protogalaxies must have something to 
do with perturbations of a typical scale-length that is larger than 
the scales responsible for the formation of primeval small pro- 
tohaloes. This lead us to the study of the influence of perturba- 
tions of various scale-lengths on the evolution of the structure. 
Preliminary results of this study very interestingly demonstrated 
that the whole supercluster- void network is caused by the inter- 
play of medium- and large-scale perturbations. These prelimi- 
nary results were reported on several conferences and summer- 
schools (see the web-site of Tartu Observatory Q). 

In the present paper, we attempt to understand the influ- 
ence of perturbations of various scale on the evolution and 
structu re o f voids in the cosmi c web. As in ISuhhonenko et al.l 
d201 ll) and lEinasto et al.1 d201 ll) . we use numerical simulations 
in boxes of various scale-lengths from 100 to 256 h~ l Mpc, cal- 
culated for power spectra cut off above different scales from 8 
to 128 h" 1 Mpc, to determine the influence of perturbations of 
various scales on the formation and internal structure of voids. 
We employ a wavelet technique to follow the evolution of under- 
dense and over-dense regions in terms of density waves of vari- 
ous scales. 

The paper is composed as follows. In the next section, we 
describe numerical models used in the study. In Section 3, we 
perform a wavelet analysis of our simulations. In Section 4, we 
describe our correlation analysis of wavelet-decomposed density 
fields. In Section 5, we investigate the structure of voids using 
haloes of various mass as objects which define voids. We discuss 
our results in Section 6, and in our last Section we present our 
conclusions. 



role. Thus, to understand the supercluster- void phenomenon cor- 
rectly, we need to perform numerical simulation in a box con- 
taining large waves. On the other hand, most systems of galaxies 
in the Universe are groups of galaxies - there are almost no very 
isolated galaxies far away from groups. The characteristic scale 
of groups is 1 h~ l Mpc, thus the simulation must have at least a 
resolution similar to this scale. 

To have both a high spatial resolution and the presence of 
density perturbations in a larger scale interval, we performed 
simulations in boxes of sizes 100 lr l Mpc, and 256 lr l Mpc, 
and resolutions A^ 3 .. = 256 3 and A^ 3 . , = 512 3 . The main param- 

gnd grid r 

eters of our series of models are given in Table Q] where L is the 
cube size, N palt is the number of particles and cells used in sim- 
ulations, and M part is the mass of a particle in units of 10 9 tf B . 
We assumed the cos mological parameters (Se liak et al.l [2005; 
iTegmarketall |2004 120061) Q m = 0.28 for the matter density, 
Qa = 0.72 for the dark energy density (in units of the critical 
cosmological density), and cr 8 = 0.84 for the initial amplitude 
parameter. We use the notation h for the present-day dimension- 
less Hubble parameter in units of 100 km s Mpc -1 ; in simula- 
tions we used a value of h = I. 

As we are interested in the study of the role of perturba- 
tions on different scales to the evolution of voids, we used sim- 
ulations with the full power spectrum, as well as with a power 
spectrum truncated at wave-numbers k cut , so that the amplitude 
of the power spectrum on large scales is zero: P(lc) - 0, if 
k < k cu t, wavelength A cut = 2jr/k cut . The cut scale in hr x Mpc, 
A cut = 2n/k cut , is given in TableQ] The amplitude of a spectrum 
was set to zero for k < k cut during the calculation of the initial 
density field, keeping all simulation parameters fixed across the 
full set of realisations. 



Table 1. Parameters of models. 



Model 


L 




AW 






(1) 


(2) 


(3) 


(4) 


M256.256 


256 


256 


256 3 


77.72 


M256.064 


256 


64 


256 3 


77.72 


M256.032 


256 


32 


256 3 


77.72 


M256.016 


256 


16 


256 3 


77.72 


M256.008 


256 


8 


256 3 


77.72 


L256.256 


256 


256 


512 3 


9.714 


L256.128 


256 


128 


512 3 


9.714 


L256.064 


256 


64 


512 3 


9.714 


L256.032 


256 


32 


512 3 


9.714 


L256.016 


256 


16 


512 3 


9.714 


L256.008 


256 


8 


512 3 


9.714 


L 100. 100 


100 


100 


512 3 


0.5583 


L100.032 


100 


32 


512 3 


0.5583 


L100.016 


100 


16 


512 3 


0.5583 


L100.008 


100 


8 


512 3 


0.5583 



2. Modelling the evolution of voids in the cosmic 
web 

Previous analyses of the observational galaxy samples and nu- 
merical simulations have shown that in the formation of super- 
clusters and voids, large-scale perturbations play an important 

1 http : //www . aai . ee/~einasto/reports . php 



Notes. 

column 1: L - size of the simulation box in h Mpc; 
column 2: A cut - cut scale in h' 1 Mpc; 
column 3: Number of particles; 
column 4: Particle mass in 10 9 M G . 

As in Suhhonenko et al. I (|20T1 . we use the notations for our 
models whereby the first characters M and L designate models 
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L256.256 
L256.128 
L256.064 
L256.032 
L256.016 
L256. 



k [h/Mpc] k [h/Mpc] 

Fig. 1. The left and right panels show the power spectra for the models of the series L256 at the epochs z — 30 and z — 0, respectively. 



with resolutions of N^a = 256 and N^d - 512, respectively. 
The subsequent number gives the size of the simulation box, 
L, in hT x Mpc, and the next indicates the maximum wavelength 
used in the simulation in h Mpc. The locations of the cells in- 
side the cubical density grid are marked by cell indices (i, j, k). 

For the models of the M256 serie s, we used in simulations 
the AMIGA code dKnebe et alJ . l200lh . This code uses an adap- 
tive mesh technique in the regions where the density exceeds 
a fixed threshold. In this code, gravity is automatically soft- 
ened adaptively, so that the softening length is near its optimum 
value in both high- and low-density regions. We chose a maxi- 
mum level of eight refinements. For models with a 512 3 resolu- 
tion, we used the GADGET-2 code with a gravitational softening 
lengt h of 10 h~ l kpc (LI 00 models) and 2 hr l kps (L256 mod- 
els) (ISpringel et all 2001: Sp ringell 120051) . The simulation LI 00 
was performed at the Leibniz-Institut fur Astrophysik Potsdam, 
and simulations M256, and L256 at the High Performance 
Computing Centre of University of Tartu. The initial density 
flu ctuation spec t rum w as generated using the COSMICS code 
by Bertschinger ( 19950; to generate the initial data, we assumed 
the baryonic matter density Qb = 0.044. Calculations started at 
an early epoch, z = 100. For every particle, we calculated the 
local density in units of the mean density, using positions of 27 
nearby particles. All models of the same series have the same 
realisation, so the role of different waves in models can be easily 
compared. To allow different models to be compared every par- 
ticle has an identification number, the same for all models of a 
series. Particle positions and velocities were extracted for seven 
epochs in the redshifts range z = 30, . . . , 0. Power spectra for the 
models of the L256 series are shown in Fig.[TJfor an early epoch, 
z = 30, and for the present epoch z — 0. 

For each particle, we also calculated the global density at the 
location of the particle (for details see Appendix A. 1 ). To find the 
global density field, we used smoothing with the B3 spline kernel 
of scale ~ 8 h~ l Mpc, which is rather close to smoothing with 
an Epanechnikov kernel of the same scale. Sm oothing of den- 
sity fie lds with different kernels is discussed bv lMartmez & Saan 
(2002). These local and global density fields were calculated for 
all models and all epochs, and were used in the subsequent ana- 
lysis to select particles belonging to a population with given 
properties, and to follow the density evolution of the model. 
Hence, for each particle we stored coordinates, local, and global 
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density values. Particles were sorted, thus their number in the 
file serves as an ID number. 

To see the effects of density waves of different scale and to 
understand the evolution of the density field, we use the wavelet 
tec hnique. We use the a t rous wavelet transform (for details, 
see lMartmez & Saarl (120021) and Appendix A. 2). The field is de- 
composed into several frequency bands as follows. The high- 
resolution (zero level) density field was calculated with the B3 
spline kernel of width equal to the size of one cell of the field, 
where every next field was calculated with twice larger kernel. 
Wavelets were found by subtracting higher level density fields 
from the previous level fields. In such a way, each wavelet band 
contains waves twice the size of the previous band, in the range 
± V2 centered on the mean (central) wave. The sum of these 
bands restores the original density field. Using this technique, 
we calculated the density fields and wavelets up to index 5 (6 for 
models of the M256 series). 

The high-resolution density fields of the model M256.256 
for the epochs z — 0, 1,5, 10 are shown in Fig. [2] The high- 
resolution density fields of models LI 00. 100 and LI 00.0 16 for 
redshifts z — and z — 2 are shown in Fig. [3] The scale of 
the cosmic web is rather different in models of different cutoff 
scale. The dependence of the scale of the web on the maximal 
wavelength of den s ity pe rturbations was investigated in detail by 
ISuhhonenko et al l (1201 ll) . 

To investigate the spatial structure of both the cosmic web 
and voids, we found haloes usi ng the adaptive Amiga Hal o 
Finder (AHF) code developed by iKnollmann & Kne be 1 2009), 
with the number of particles in a halo A^ p > 20. Haloes and their 
parameters (masses, virial radii, positions, velocities etc.) were 
found for all models and simulation epochs. 

3. Analysis of models 

To understand the evolution of the density field, and to see the 
effects of density waves of different scale, we analyse the den- 
sity fields and their wavelet decompositions for various simula- 
tion epochs and cut off scales. We focus on the evolution of the 
following properties: 

1 . the global patterns of the density field, using wavelets of dif- 
ferent levels; 

2. the fine structure of the density field for perturbations of var- 
ious scales; 

3. density perturbations of various scales; 
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Fig. 2. The high-resolution density fields of the model M256.256 are shown in the left column, at the k = 140 coordinate. The 
second, third, and fourth columns show the wavelet w5, w4, and w3 decompositions at the same k, respectively. The upper row gives 
data for the present epoch, z — 0, the second row for the redshift z — 1, the third row for the redshift z — 5, and the last row for 
the redshift z = 10. Densities are expressed on a linear scale. In wavelet panels, green and red colours show the positive regions of 
wavelets, and the blue colour shows negative (under-density) wavelet regions. Colour codes are plotted at the bottom of the figure 
for wavelet w3 at epoch z — 10. 



4. the density distribution in void regions. 

3. 1 . The evolution of the global patterns of density fields 

To compare the patterns of density fields and their wavelet 
decompositions for various cosmic epochs, we use the model 
M256.256. 

In Fig. [2] we plot the high-resolution density field and 
wavelets of the model M256.256 at the four redshifts: z = 
0, 1, 5, 10. In the second column of Fig. [2] the wavelets of 
order five are shown. The characteristic scale of density pertur- 



bations for this wavelet is 64 h~ x Mpc. The upper density levels 
used in plotting for w5 are 1.2, 0.7, 0.2, and 0.1, for redshifts 0, 
1, 5, and 10, respectively. Lower limits to these redshifts are 
-0.6, -0.35, -0.1, and -0.05. In wavelets of the order4 and 3, a 
similar choice of colour limits is applied. This colour-coding of 
wavelets at different redshifts is chosen so that a certain colour 
corresponds approximately to the density level, corrected by the 
linear growth factor for that redshift. Wavelet blue colours cor- 
respond to under-dense regions of density waves, green colours 
to slightly over-density regions, and red colours to highly over- 
density regions. 
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Fig. 3. Zoom-ins to the high-resolution density fields of the models L100.100 and L100.016, left and right columns, respectively. 
Zoom factor is 2, central 50 x 50 hr x Mpc (256 x 256 pixels) of all models are shown. Upper panels are for the present epoch 
z — 0, lower panels for the epoch z - 2. All panels are at the k = 51 coordinate. Cross sections (beams) at coordinates j = 222 
and k = 51 for both models are shown in Fig. [5] for three redshifts to see the evolution of the density field and its wavelets. In the 
upper left corner of the Figure there is a rich supercluster i n the model L100. 100, absent in the model L100.016. Both models have 
at the right edge of the Figure a rich cluster. This cluster is well seen in the Fig. [5] at i = 390 coordinate. Densities are expressed 
in the logarithmic scale, identical lower and upper limits for plotting with the SAO DS9 package are used. The border between the 
light blue and the dark green colo urs corresponds to the critical dens ity Di oc = 1 .6, which separates low-density haloes and haloes 
collapsed during the Hubble time (Kaiser, 1984; Bard een et allll986h . Note that in both models and simulation epochs the majority 
of filaments in voids have densities below the critical density. 



Figure|2]shows that the pattern of the cosmic web on wavelet 
w5 is almost identical at all redshifts, only the amplitude of the 
density waves increasing approximately in proportion to the lin- 
ear growth factor. This linear growth is expected for density 
waves of large scales, which are in the linear stage of growth. 
The pattern of the web of the wavelet w4 changes little, but the 



growth of the amplitude of density waves is more rapid. The pat- 
tern of the wavelet w3 changes much more during the evolution, 
and the amplitude of density waves increases more rapidly, but 
essential features remain unchanged, i.e. the locations of high- 
density peaks and low-density depressions are almost indepen- 
dent of the epoch. 
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Fig. 4. Wavelets wl of models L100.100 and L100.008 at redshift z — 10 are shown in the left and right panels, respectively, at 
coordinate k — 51. Densities are expressed on a linear scale, only over-density regions are shown. As in Fig. [3] the central 256 x 256 
pixels of the full models are shown. The characteristic scale of density perturbations corresponding to this wavelet is 0.4 hr x Mpc. 
Note the weakening of peak densities of the model L100.100 in the region of future large under-dense region, seen in Fig. |3] 



Fig. |2] shows that at all redshifts high-density peaks of 
wavelets of medium and large scales almost coincide. In other 
words, density perturbations of medium and large scales have a 
tendency of phase co upling or synchronisation at peak positions. 
Eina sto et al.l d201 lb reached the same conclusion using models 
with a much broader scale interval. Figure |2]shows that the syn- 
chronisation of medium and large scales applies also to under- 
dense regions. The analysis below describes the differences be- 
tween the synchronisation of over- and under-dense regions in 
quantitative terms. 



3.2. The evolution of the fine structure in the density field 

To follow the evolution of the fine structure in the density field, 
we compare high-resolution density fields of the models of LI 00 
series (see Fig. [3}. We show slices at the coordinate k — 51 of the 
full model L100. 100, and of the strongly cut model L100.016, at 
epochs z — 0, and z — 2. The k coordinate is chosen so that 
the slice in the model L100.100 crosses a large under-dense re- 
gion between a rich supercluster and several rich clusters. To see 
more clearly the differences between the density fields of models 
L100.100 and L100.016, we show in Fig.[3]only the zoom-in of 
the central 50 x 50 h Mpc (256 x 256 pixels) region. To com- 
pare the present field with the initial density field, we use the 
smallest scale wavelet wl for both models at the redshift z = 10, 
shown as a zoom-in plot in Fig.|4] which is similar to the plot for 
the present epoch in Fig. [3] 

The power spectrum of density perturbations has the highest 
power on small scales. Thus, the influence of small-scale per- 
turbations relative to large-scale perturbations is strongest in the 
early period of evolution of structure. For this reason, density 
fields and wavelets wl at early epochs are qualitatively rather 



similar for the full model LI 00. 100, and for the model cut on 
small scales, A cM = 8 h~ x Mpc, called L100.008 (see Fig.g}. 

However, there are small but important differences in the 
patterns of small-scale structures in the models LI 00. 100 and 
L100.008 at z = 10. The density peaks of the model L100.008 
have more or less equal heights throughout the whole simulation 
box, whereas in the model LI 00. 100 in regions of future voids 
the peak heights are lower than in the future supercluster regions. 
This shows that already at the epoch z — 10 large-scale pertur- 
bations started to influence the density field on small scales. 

The colour coding is identical in all panels of Figs. [3] Codes 
are chosen so that the density level Di oc = 1.6 is clearly visible. 
Particle clouds of density above this limit can form collapsed 
haloes during the Hubble time, and below this limit they stay in 
pre-g alactic more diffuse form. The act ual collapse parameter is 
1.69 dKaiseAfl984T:lBardeen et aUll986l) . In both models, there 
exists a network of filaments. We see that in all models shown in 
Fig. [3] the majority of filaments have densities below this limit, 
i.e. that filaments consist of strings of primordial matter that have 
not yet formed compact haloes. 

At the epoch z — 2, there are already strong differences be- 
tween the density fields of models L100.100 and L100.016. In 
the model L100.016, the high-density knots are distributed more 
uniformly, and the whole pattern of filaments has a smaller scale. 
We discuss the distribution of mean void sizes for both mod- 
els below. In the full model L100.100, the fraction of particles 
above the critical density D] oc = 1 .6 in regions of low global 
density (voids) is lower (see Fig. |7). Regions containing high- 
density knots (marked by green and red) start to form superclus- 
ters. In the model L100.016, these regions are more uniformly 
distributed. 

The comparison of distributions for the epochs z — 2 and 
z — shows that filaments contract with time, as known from 



7 



J. Einasto et al.: Void structure in the cosmic web 



L 1 00. 1 00 222,5 1 z=0 LI 00.0 16 222,51 z=0 




i [0.195»Mpc/h] i [0.195*Mpc/h] 



Fig. 5. The evolution of the local density and wavelets of the models L100.100 (left panels) and L100.016 (right panels) in beams 
along the /-coordinate at j = 220, k — 51. The same ^-coordinate was used in plotting the density field in Fig. [3] Data are shown 
for epochs z — 0, 2, 5. To see better details only the region 144 < i < 400 of length 50 h~ l Mpc is shown. The characteristic scale 
of the wavelet w5 is 12.5 hr l Mpc. Wavelets are divided by the factor / oc (1 + z) _1 . 



many earlier studies cited above. In the model L100.016, the 
pattern of filaments changes very little between the epochs z — 2 
and z — 0. In contrast, for the model L100. 100 density evolution 
can be clearly seen. At z = 2, small-scale filaments fill almost the 
whole space of voids between rich systems. At z — 0, most of 



these small filaments have merged leaving more space for very 
low-density regions. 

The differences between the models L100.100 and L100.016 
at the present epoch z = are very well seen in Fig. [3] In the 
model LI 00. 100, between the supercluster at the left corner and 
the cluster at the right edge there is a large low-density region. 
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Fig. 6. The local density and wavelets of the models L256.256 (left panels) and L256.008 (right panels) in beams along the 
/-coordinate at j = 90, k = 105 at the early epoch, z = 30. To enhance details, only the region 100 < i < 250 of length 
75 hT l Mpc is shown. The characteristic scale of the wavelet w5 is 32 hr x Mpc. 



This region is crossed near the center by a filament that has sev- 
eral knots in the green and red colour. In the model LI 00.0 16, 
there are no rich superclusters, the whole region being covered 
by a web of small-scale filaments. In other words, large-scale 
perturbations present in the model L100.100, have suppressed 
the growth of the density of filaments in void regions. 

There exists a low-density smooth background, seen in the 
Fig.[3]in deep dark-blue colour. The density of this background, 
D * 0. 1 , is lower at the present epoch z — 0, i.e. the density of the 
smooth background decreases with time (see Fig. |7j. Regions of 
very low density have much larger sizes in the model L100.100 
than in the model L100.016. 

3.3. The evolution of density perturbations of various scales 

We now follow the evolution of density perturbations of vari- 
ous scales in the models of the series LI 00 and L256 in more 
detail. The evolution of density perturbations has several char- 
acteristics: the shapes of density perturbations of various scales 
and their change with time; synchronisation, amplification, and 
suppression of density perturbations of various scales; and the 
formation of regions of very low density. 

To see the evolution of the density field and its wavelets, 
we show in Fig.|5]one-dimensional density distributions (beams) 
along horizontal lines of the plane in Fig. [3] The distributions 
are generated along the axis j = 222, k = 51, for the models 
L100. 100 and L100.016 at the redshifts 0, 2, 5. Beams are taken 
along the /-coordinate; at each i value, all cells in the j- and 
^-coordinate within ±5 from the centre of the beam are counted, 
i.e. beams have sizes 11x11 cells (about 2.15X2.15/T 1 Mpc 
in the models of the L100 series). The average densities for a 
given i are found, and are given in the mean density units. We 
use overdensities D — 1 here, so the mean overdensity level is 
zero, which is similar to the mean value for wavelets. Wavelet 
amplitudes are divided by the linear growth factor, thus during 
linear growth their amplitudes should not change. The effective 
scale for the wavelet w5 of the model L100 is 12.5 h' 1 Mpc, as 
seen also from the separation between the high-density peaks in 
Figs.|5] 

At high redshifts, perturbations of various scales have almost 
sinusoidal shapes, and the wave peaks have approximately equal 



heigh ts. During the subsequent stages of evolution. lEinasto et al.l 
(1201 ll) showed that perturbations of larger scales begin to affect 
the evolution. These perturbations amplify small-scale perturba- 
tions near maxima, and suppress small-scale perturbations near 
minima. Thereafter, still larger perturbations amplify smaller 
perturbations near their maxima, and suppress smaller pertur- 
bations near their minima, and so on. 

For early stages, the density fields and wavelets are given for 
the models L256.256 and L256.008, for the epoch z = 30, in 
Fig. [6] This model also allows us to illustrate the influence of 
larger waves, because the characteristic scale for the wavelet w5 
of the model L256 is 32 h~ x Mpc. 

We see that in the model L256.008 wavelets of different or- 
der have approximately sinusoidal shapes. The wavelet w5 of 
this model has zero amplitude, and the amplitude of the wavelet 
w4 is small, because in this model only waves of scale up to 
8 /i~'Mpc are present. The wavelet w3 of the effective scale, 
which is approximately equal to the cut scale of this model, has 
the highest amplitude. The maxima and minima of this wavelet 
are partly synchronised with the maxima and minima of the 
wavelets w2 and wl; synchronisation is better in the regions that 
happen to coincide with the maxima and minima of the wavelet 
w4. 

In the model L256.256 at the early epoch z = 30, the wavelet 
w5 has the highest amplitude, the next highest being that of the 
wavelet w4. These two wavelets are only partly synchronised, 
i.e. the maxima of the wavelet w5 coincide in position only ap- 
proximately with the maxima of the wavelet w4, and near the 
minima of the wavelet w5 there is a secondary maximum of the 
wavelet w4, i.e. the wavelet w4 behaves as the first overtone of 
the wavelet w5. The shape of small-scale wavelets is not sinu- 
soidal. This is caused by large-scale waves that have started to 
change the shapes of small-scale waves. 

Further evolution with time of the wavelets can be followed 
using the model L100 in Fig. We discuss first the evolution 
of wavelets in the model L100.016 (see the right panels of the 
Figure). The largest wavelet w5 has for all redshifts approxi- 
mately the same amplitude and a sinusoidal shape, suggesting 
that density perturbations of this scale are in the linear growth 
regime. The shape of the next wavelet w4 is very different from 
sinusoidal. In some regions of maxima of the wavelet w5, the 
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wavelet w4 has very strong maxima. An example is the region 
near i ss 370, where all wavelets of smaller scale also have 
strong maxima, and wavelets of all scales up to w4 are very 
well-synchronised. The overall shape of the density profile is 
determined by the wavelet w3 that has a characteristic scale 
3.1 h~ l Mpc. Most peaks seen in the density profile are due to 
the maxima of this wavelet. In most of these peaks, wavelets of 
smaller scale also have maxima, i.e. near the peaks small-scale 
wavelets are synchronised. 

When we compare the density and wavelet distributions of 
this model at various epochs, we see little difference. This shows 
that in the model L100.016 the structure has rapidly evolved at 
early epochs and changes only a little later. The most important 
development is the decrease in the density in deep void regions, 
and the increase in the density in massive haloes. 

We now consider the evolution of the full model L100.100, 
shown in the left columns of Fig. [5] We see that the evolution 
of the largest wavelet w5 is almost linear up to the epoch z > 1 
(the density and wavelet distributions for the epochs z = 1 and 
z — 2 are rather similar, only the amplitudes of wavelets up to w4 
are smaller for z = 1). The shapes of the wavelet w5 for differ- 
ent epochs are almost sinusoidal and the heights of the maxima 
are approximately equal. The next wavelet w4 behaves as a first 
overtone of the wavelet w5 - near the minima of w5 there are 
maxima of w4, which have much lower heights than the maxima 
near the maxima of w5. Evidence of this phenomenon can also 
be clearly seen in the wavelet analysis of the Sloan Digital Sky 
Survey (see Fig. 6 of Einasto ^t al.l (1201 lh ~). Near the joint max- 
ima of w5 and w4, there are very strong maxima of all wavelets 
of smaller order; this is very well seen at locations i w 290 and 
i ~ 380. 

Near the minima of w5 and the maxima of w4, wavelets of 
smaller order also have maxima, but these maxima get weaker 
at lower redshifts, as in the region around i « 220. Here, small- 
scale wavelets are partly synchronised, and small-scale peaks of 
the density field are related to maxima of the wavelets w2 and 
w3. 

The density and wavelet distributions of the model LI 00. 100 
for the present epoch z = are completely different from the dis- 
tributions at higher redshifts. In the present epoch the dominant 
feature is a large under-dense region in the interval 120 < i < 
380. This large under-dense region is caused by large-scale den- 



sity perturbations that are not shown as wavelets in the Figure. 
In this region, these large-scale density waves have their min- 
ima and the amplitudes of density waves of smaller scales, in- 
cluding w4 and w5, are suppressed. The amplitudes of wavelets 
w3 and lower orders are almost zero. Near the density maxima 
seen at higher redshifts at i ~ 210 and i m 290, there are very 
weak density peaks with maxima below the mean density level. 
These maxima are seen in the density field as weak filaments 
(see Fig. [3]). 

The most remarkable feature of the density field in the model 
L100.100 at the present epoch is the presence of large under- 
dense regions of very low density D « 0.1, which can be clearly 
seen in Fig. [3] in deep-blue colour. At earlier epochs, the den- 
sity in these regions was higher and there were numerous low- 
density peaks within the regions, for the present epoch most of 
these peaks are gone. This is caused by density perturbations on 
larger scales. 

When we compare the evolution of density distributions 
of models L100.100 and L100.016, we see remarkable differ- 
ences. These differences are solely due to the presence of den- 
sity perturbations of scales larger than A cut = 16 h~ l Mpc in the 
model LI 00. 100. We note that both models were generated with 
identical "random amplitudes", i.e. the perturbations of scales 
A < Awn = 16 h Mpc are identical in both models. 

The main result of the paper by Eina sto et al.l (1201 ll) was 
that at all redshifts high-density peaks of wavelets of large and 
medium scales almost coincide. Figures [2] [3] and [5] show that 
the same conclusion is valid for positions of density depres- 
sions (deepest voids) of wav elets of large and m edium scales. 
The other main conclusion of lEinasto et al.l (1201 ll) was that po- 
sitions of peaks of waves of different scale coincide, i.e. density 
waves of different scale are synchronised. A look at Figs. |2]and 
[5] shows that the synchronisation of density waves of different 
scales also concerns density depressions of large and medium 
scales. However, the synchronisation of the depressions of den- 
sity waves, which is responsible for the formation of voids, is 
less pronounced than that of density peaks. 

3.4. The evolution of density distributions in void regions 

To investigate the evolution of the density distribution in void 
regions, we selected in all models void particles in the present 
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epoch z — 0. We calculated the distributions of the global densi- 
ties of full models (i.e. models with no cuts in the power spec- 
tra), and selected particles of the lowest global density values, 
about 10 % of all particles (13.390.895 particles in the model 
L100. 100). The corresponding value of the threshold global den- 
sity is 0.565 for the model L100.100 (in units of the mean den- 
sity). 

The evolution of the distributions of local densities of void 
particles for the full model LI 00. 100 is given in the left panel of 
Fig. [7] The distributions of void particle local densities for the 
models of the LI 00 series for the present epoch are shown in the 
right panel of Fig. [7] 

Figure [7] shows that the initial distribution of particles at the 
redshift z - 10 is quite symmetrical in the log-log representation. 
The density distribution has a peak at D « 0.8. As time goes on, 
the peak density decreases, and at the present epoch has a value 
D as 0.2; the lowest density occurs close to D m 0.1. Particles in 
overdense regions (D > 1) form haloes; with decreasing redshift 
z, the peak densities increase, i.e. haloes in void regions become 
more massive and denser. 

Density evolution depends strongly on the cutoff wavelength 
of the model. In the most strongly cut model L100.008 at the 
present epoch, the fraction of particles in very low-density re- 
gions (D at 0.2) is about five times lower than in the full model 
L100.100. In the model L100.008, most particles in void regions 
form haloes of mean density D as 100, and with the maximum 
density D as 1000. With the increase in the cut wavelength A cut , 
the fraction of particles in very low-density regions increases. 
The maximal density of haloes in void regions reaches the high- 
est value, D as 2500, in the model L100.016. In this model, 
density waves between the scales 8 and 16 /i 1 Mpc amplify the 
density of haloes. If density perturbations of larger scale are in- 
cluded (models L100.032 and L100.100), then in the void re- 
gions still larger perturbations start to decrease the maximum 
masses of haloes. This depression is largest in the full model 
LI 00. 100, where the maximum local densities in void haloes 
reach the values D « 600. 



4. Correlation analysis of wavelet-decomposed 
density fields 

We now attempt to quantify some of the qualitative statements 
given earlier in the text. To this end, we perform the correlation 
analysis of wavelet-decomposed density fields. Our a pproach is 
analogous to the one presented in Ein asto et al.l (1201 ll) with the 
exception that here instead of over-densities we focus on under- 
dense regions. In the following, we present our results only for 
the model M256 since the other ones lead to very similar results. 
We consider six wavelet levels: wl, w2, . . ., w6 with the effective 
smoothing scales of 4, 8, . . ., 128 h Mpc, respectively, and use 
the simulated density fields at five different redshifts z = 30, 10, 
5,1,0. 

Quite generally, one can choose two redshifts, z, and Zj, 
along with two wavelet levels, w,„ and w n , and calculate the cor- 
relators 

{(6 WmZi - (6 WrnZi ))(6 w „ Zj - (6 KiiZj ))) 

t~W m Zi,W„Zj i ' \i) 

y<(<W f - <<W) 2 X(<Wj - <<W) 2 > 

where 6 WmZi corresponds to the wavelet-decomposed density field 
for level w m at redshift z,. The angle brackets represent an en- 
semble average, which under the ergodicity assumption is re- 
placed by a simple spatial average. We note that we calculate the 



correlators for zero lag only, i.e. we do not shift one field with 
respect to the other. 




wavelet level 



Fig. 8. The behaviour of the correlation coefficient for different 
redshift pairs ( Zi ,Zj) e {(30, 10); (30, 5); (30, 1); (30,0)} for all 
the six wavelet levels. The solid and dashed lines correspond 
to the under- and overdensities, respectively (10% of the most 
under/over-dense cells). We see that for the largest smoothing 
scale, i.e. w6, all the correlators stay quite close to r — 1, while 
later on, as the other redshift decreases below z = 30, the lines 
start to deviate from r = 1 . Thus, on the largest scales the infor- 
mation is approximately preserved, while on the smallest scales 
the information gets gradually erased. 

In what follows, the above correlators are not calculated for 
the full density fields, but instead masks are applied to separate 
under-dense void regions. If two wavelet fields S WmZj and 5 WnZ . 
are correlated then the mask is always defined with the field that 
has larger smoothing scale, e.g., the field S WmZl if w m > w„. The 
masking level is taken such that only 10% of the most under- 
dense cells end up inside the mask. 

In the following, we use two types of correlators: 

1. Fixed wavelet scale correlators, i.e. w m = w n , at different 
redshifts. 

2. Correlators at fixed redshifts, i.e. z, = Zj, for different wavelet 
levels. 

In the first case, we take one of the density fields always at 
redshift z = 30, which is high enough for all of the scales of in- 
terest to be well in the linear regime. It is easy to understand that 
under the linear evolution, where the values of 6 just get multi- 
plied by the same scale-independent but time-varying factor, the 
correlation coefficient should always stay at the value r = 1, i.e. 
all of the initial information is well preserved. In Fig. [8] we show 
the behaviour of the correlation coefficient for different redshift 
pairs: ( Zi ,Zj) e {(30, 10); (30, 5); (30, 1); (30,0)} for all the six 
wavelet levels. It is easy to see that for the largest smoothing 
scale, i.e. w6, all the correlators stay rather close to r = 1, while 
later on, as the other redshift gets smaller than z = 30 the lines 
start to decline from r = 1, especially on the smallest scales. 
Thus, on the largest scales the information is approximately con- 
served, while on the smallest scales the information gets erased. 
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wavelet level 

Fig. 9. Correlators at fixed redshifts Zi =Zj = 30, 10, 5, 1 , from 
top to bottom (from light to dark blue) for the under-dense re- 
gions. For w m = 1, the curves are peaked at w„ = 1, while they 
drop gradually as w n is increased to higher values. Similarly, the 
curves for w m = 2 are peaked at w„ = 2, and get reduced as the 
distance increases from this point. For the other w m and w„ val- 
ues, the behaviour is very similar. We see that the lower the out- 
put redshift, the narrower the coupling kernels, i.e. a nonlinear 
evolution in under-dense regions leads to additional decoupling 
of nearby wavelet modes. 



The lower the redshift of the other density field the greater the 
loss of information. In practice, for the cases z — 10 and z = 5 
the loss of information is relatively modest if the wavelet level 
w > 3. For z = 1 and z — 0, the information is approximately 
saved only for the largest scales. The dashed lines in Fig.|8]show 
the corresponding results for the over-dense regions, in this case 
focusing on the 10% of the most over-dense cells. As one can 
see, the information loss in under-dense regions occurs more 
rapidly with time than for the over-densities. 

In Fig. [9] we plot the correlators at fixed redshifts zi = Zj - 
30, 10, 5,1,0 from top to bottom (from light to dark blue). For 
m = 1, the curves are peaked at n = 1, while their amplitude 
decrease gradually as n increases to higher values. Similarly, the 
curves for m — 2 are peaked at n = 2, and decrease in ampli- 
tude as one moves to neighbouring wavelet levels. For the other 
n values, the behaviour is very similar. As long as the evolu- 
tion proceeds in a linear manner, i.e. the growth depends only 
on redshift, but is independent of the wavelet scale, the coupling 
kernels plotted in Fig.|9]should stay exactly the same. However, 
we see that the lower the output redshift, the narrower the cou- 
pling kernels, i.e. a nonlinear evolution in under-dense regions 
leads to the additional decoupling of the nearby wavelet modes. 

The correspon ding figure for the over-densities was given in 
lEinasto et al.l (1201 ll) (see Fig . 7 the re). The main difference be- 
tween Fig. 7 bv lEinasto et al.l (1201 ll) and Fig. [9]is that in the case 
of over-densities the coupling kernels become broader because 
of nonlinear evolution, i.e. instead of desynchronisation we have 
increasing synchronisation of over-densities of different wavelet 
levels. 

It is important to realise that even with only the linear evo- 
lution of the Gaussian density field the nearby wavelet levels 
at fixed redshift get significantly coupled, since the neighbour- 
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Fig. 10. The cumulative mass functions of the AHF haloes of 
models of the L100 series with various cutoff scales. The up- 
per, middle, and lower panels are for the redshifts z — 0, 2, 5, 
respectively. 



ing levels tend to contain some of the common Fourier space 
modes. However, assuming only linear evolution it is clear that 
the coupling does not change with redshift. 



5. The structure of voids 

Cosmic voids are defined by objects surrounding them - galax- 
ies and clusters of galaxies of various luminosity (mass). In the 
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Fig. 11. The left panels show the numbers of voids, defined by the AHF haloes for various threshold masses and models, as shown 
in Fig. [10] The right panels show the mean radii of voids, defined by the AHF haloes for different threshold masses, M t h, and for 
various cut-off scales. The upper, middle, and lower panels are for the redshifts z = 0, 2, 5, respectively. 



case of models, it is customary to use dark matter haloes instead 
of galaxies or clusters. To investigate the influence of density 
perturbations of various scale on the void structure, we shall use 
our models of the L100 series which have the highest resolution 
in mass a nd scale. To find haloes, w e applied the Amiga halo 
finder by iKnollmann & Knebd d200S3t) . We characterise the web 
and void structure by the dependence of the halo mass function 
on the scale of density perturbations for various epochs, and by 
the number and radii of voids defined by haloes of various mass. 
Our results are shown in Figs. HTJIandHTI 



Fig. [10] shows the cumulative mass functions of the AHF 
haloes for all models of the L100 series, for three epochs, z = 
0, 2, 5. We see that in the models where the large-scale waves 
have been cut off, maximum haloes have much lower masses 
than in the full models. This effect can also be seen in Fig. [3] 
The differences between the models with various cutoff scales 
increase with time: at early epochs, halo masses are lower and 
time is needed for most massive haloes to grow. 

To find voids, w e used a simple void finder proposed by 
Einas to et a i1 (fl989h . For each vertex of the simulation grid, we 
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first calculated its distance to the nearest AHF halo. The max- 
ima of the void distance matrix correspond to the centres of 
voids, and their values are the void radii. The distribution of the 
AHF haloes is noisy, thus there are many nearby local maxima 
in the distance matrix. We define the position of the void centre 
as the location of the cell, which has the largest distance to a halo 
within a box of the size of +3 grid elements. 

The left panels in Fig. Qj] show the number of voids found 
for various AHF halo mass thresholds, which correspond to sys- 
tems of galaxies of different mass. The right panels of the Figure 
show mean radii of voids as a function of the AHF halo mass 
threshold. The void numbers and radii characterise the hierar- 
chy of voids. We see that as the threshold halo mass used in 
the void search increases, the number of voids continuously de- 
creases, and the void radii increase. This means that some fila- 
ments are fainter than the respective mass threshold limit, and 
do not contribute to the void definition. Both parameters also de- 
pend strongly on the highest density perturbations used in the 
models. Models cut at larger scales have more voids, but their 
radii are smaller. This effect can be clearly seen at all simulation 
epochs. 

In the models LI 00. 100 and LI 00.032, the dependence of 
the number of voids and their mean radii is a more-or-less con- 
tinuous function of the halo threshold mass. In contrast, in the 
strongly cut models L100.016 and L100.008, at higher mass 
thresholds, the number of voids decreases very rapidly and void 
radii also increase more rapidly than in the models LI 00. 100 and 
L100.032. This effect is due to the very sharp decrease in the 
number of haloes of high mass (see Fig. ITOb . These rare haloes 
define very large voids. The sizes of these voids are not charac- 
teristic of the overall cosmic web pattern of the particular model. 

The simulation box used in the present void structure study 
has the size L - 100 lr l Mpc, thus very large perturbations 
responsible for the formation of rich superclusters are absent 
in the model. For this reason, the largest voids for the highest 
halo mass thresho l ds ha ve radii 10 hr l Mpc. As shown by 
Suhh onenko et al.l (1201 ll) . in models with cube sizes L = 265 
and L = 768 h~ l Mpc the maximum v oid radii are much larger 
(see Fig. 6 by Suhh onenko et al.l (120 1 ll) ). This comparison shows 
that large-scale density perturbations are needed to form voids 
defined by superclusters. 

6. Discussion 

According to the current cosmological paradigm, the cosmic 
web with systems of galaxies of various scale and mass, from 
clusters to filaments and superclusters, and voids between them, 
is formed from tiny density perturbations during the very early 
stage of the evolution of the Universe. For the formation of the 
web and of voids between various objects of the web, the pres- 
ence of a continuous spectrum of density perturbations of vari- 
ous scales is essential. The power spectrum of density perturba- 
tions has the highest power on small scales. Thus, the influence 
of small-scale perturbations relative to large-scale perturbations 
is strongest in the early period of structure evolution. Small-scale 
systems, i.e. small haloes in simulations and dwarf galaxies in 
th e real world are the earlies t compact objects to form. As shown 
bv lSuhhonenko et al.l (120 1 ll) . and confirmed in the present study, 
small-scale halo es form at early epoch s everywhere. The wavelet 
analysis done by Einast o et al.l (1201 ll) shows that wavelets wl at 
redshift z — 30 are almost identical in models L256.256 and 
L256.008 (see Fig. 8 there). The present study shows that this 
is also the case for the wavelets wl of the models L100.100 and 
L100.008 at the epoch z = 10 (see Fig.g}. 



Further evolution of the web depends on the presence of den- 
sity perturbations on larger scales. In models where medium- 
scale density pe rturbations are absent, no systems of filaments 
and voids form ( Suhhonenko et all 1201 ll) . In models with den- 
sity perturbation spectra cut on large scales, the cosmic web with 
filaments and voids has a characteristic scale of the largest scale 
present in the density perturbation field. The main quantitative 
characteristics of the web - the masses of haloes, the number 
and sizes of voids defined by haloes of various mass - depend 
on the largest scale perturbations present. 

In the models with strongly cut power spectra (A cut < 
16 h' 1 Mpc), the maximum masses of haloes are lower than in 
the models with a full power spectrum. Their number is larger, 
but they define a cellular cosmic web with smaller mean void 
sizes (see Fig. [TTJ. With the increase of the cut-off wavelength 
/l cut , the maximum masses of haloes increase, and they define a 
cellular web with larger cells but fewer voids. 

The wavelet analysis described in previous Sections shows a 
very important property of the evolution of density waves with 
time: the sy nchronisation of the p hases of density waves on var- 
ious scales. lEinasto et al.l (1201 ll) discussed this property of the 
evolution of over-density features - clusters and superclusters of 
galaxies in numerical simulations. In the present paper, we have 
followed the evolution of both over- and under-density regions. 
The analogy in the evolution of over- and under-density regions 
is expected, since in the early linear stage of the evolution of 
structure positive and negative parts of density waves were sim- 
ilar and symmetrical. 

The wavelet analysis leads us to the conclusion that the prop- 
erties of the large-scale cosmic web with filaments and voids 
depend on two connected properties of the evolution of density 
perturbations. The first property is the synchronisation of density 
waves of medium and large scales. Due to the synchronisation of 
density waves of different scales, positive amplitude regions of 
density waves add together to form rich systems of galaxies, and 
negative amplitude regions of density waves add together to de- 
crease the mean overall density in voids. The amplification of 
density perturbations is another property of density evolution. 
Due to the addition of negative amplitudes of medium and large 
scale perturbations, there is no possibility for the growth of the 
initial small-scale positive density peaks in void regions. For this 
reason, small-scale protohaloes dissolve there. In the absence 
of medium and large-scale density perturbations, these peaks 
would contract to form haloes, which would also fill the void 
regions, i.e. there would be no void phenomenon as observed. 

Simulat i ons w ith tr uncated power spectra were performed by 
iLittle et all (119911) and lEinasto & Gramannl (1 19931) . iLittle et al.l 
(1991) used three-dimensional simulations of resolution 128 3 
with power spectra in the form P(k) ~ k~ l for k < k c and 
P(k) = for k > k c , i.e. the spectra were cut at small scales, 
in contrast to our study here. The cuts were scaled so that 
k — 1 represents the fundamental mode of the simulation box 
of the size L — 64 ft -1 Mpc. The authors used the cuts: k = 
2, 4, 8, 16, 32, 64. The main result of the study was that the 
structure of the cosmic web depends on density perturbations of 
larger scale than the cut-off scale, in accordance with our results. 
lEinasto & Gramannl d 19931) made a two-dimensional simulation 
of the resolution 51 2 2 with the full power spectrum, and a power 
spectrum cut at the scale A, = L/4, where L is the size of the 
simulation box. The fine structure of filaments in both models 
was rather similar, only the location and strength of filaments 
was slightly different. This result is also in agreement with our 
present findings. 
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As the analysis shows, the phase synchronisation of both 
positive and negative sections of density waves is stronger 
for density waves of larger scales, A > 32 h~ x Mpc. Scales 
larger than the sound horizon at recombination, m 146 Mpc ac- 
cordin g to the most recent cosmological data by Jarosik et al.l 
(2010), were outside the horizon most of the time. This scale, 
105 hT x Mpc for the presently accepted Hubble constant h = 
0.72, is surprisingly close to the charac t eristic scale of the 
supercluster-void network (Eina sto et al.l Il997al 1200 lb . The 
skeleton of the supercluster-void network was created during the 
very early post-inflation stag e of the evolution of the Universe 
dKofman & Shandarinl Q988). This result is also true for large 
voids between superclusters of galaxies - the seeds for these su- 
pervoids were created in the very early Universe. 

7. Conclusions 

Our present study of the evolution of density perturbations of 
various scales has led to the following conclusions: 

- The formation of the cosmic web with filaments and voids is 
due to the synchronisation of density waves of medium and 
large scales, and the amplification of both over- and under- 
dense regions. 

- Voids are regions in space where medium- and large-scale 
density waves combine in similar under- density phases. 

- Owing to phase synchronisation, the mean density of matter 
in void regions is below the mean density, thus initial small- 
scale perturbations cannot grow. 
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Appendix A: Density field and wavelets 

A.1. Density field 

For each particle, we calculated the global density at the loca- 
tion of the particle. For this purpose, we first found the high- 
resolution density field, using a B3 spline 

B 3 (x) = T- [l* ~ 2| 3 - 4|x - 1| 3 + 6\x\ 3 - 4\x + 1| 3 + \x + 2| 3 ] ; (A. 

this function differs from zero only in the interval x e [-2,2]. 
The one-dimensional B 3 box spline kernel of width h — N is 

K$\r, N) = B 3 (x/N)/N. (A.2) 

This kernel preserves the interpolation property (mass conserva- 
tion) for all kernel widths that are integer multiples of the grid 



step, h = N. The 3-D K B box spline kernel we use is given by 
the direct product of the three one-dimensional kernels 

K B (x; N) = K (3 \x; N) = K^)(x; N)K^\y; N)K { *\z; N), (A3) 

where x = {x, y, z}. To calculate the high-resolution density field, 
we use the kernel of scale, equal to the cell size of the particular 
simulation. 

A.2. Wavelets 

We use the a trous wavele t transf orm (for details see lStarck et all 
I1998I; IStarck & Mur tagh, 2002]). The field is decomposed into 
several frequency bands as follows. The high-resolution (zero 
level) density field was calculated with the B 3 spline kernel with 
width equal to the size of one cell of the field, every subsequent 
field being calculated with a kernel twice as wide. Wavelets were 
found by subtracting higher level density fields from the previ- 
ous level fields. In such a way, each wavelet band contains waves 
twice the scale of the previous band, in the range ± V2 centered 
on the mean (central) scale. The sum of these bands restores the 
original density field. 

The 'a trous algorithm' wavelet transform decomposes an 
nxnxn data set D as a superposition of the form 

,/ 

D = Dj + J] wj, (A.4) 

7=1 

where Dj is a J times smoothed version of the original data D, 
and wj represents the structure of D at scale 2 J . The wavelet 
decomposition output consists of / three-dimensional mother 
fields Dj and wavelets Wj of size n x n x n. Following the tra- 
ditional indexing convention, we mark the mother fields and 
wavelets of the finest scale with the index j = 1. The smoothed 
version of the original data, Dj — Do, is the density field found 
with the kernel of the scale, equal to the cell size of the simula- 
tion L/Vgrid. 

The wavelets can be found in a recursive manner, but we 
also needed to evaluate the partial density fields (mother fields 
of different order) for our analysis. Thus, we found the mother 
fields Dj first by convolving the field by the B 3 kernel of 
twice the scale used for calculating the field . We then found 
the wavelets of index j by subtracting the mother density fields 

Wj = Dj_i - Dj. (A.5) 

In this construction, a wavelet of index j describes density waves 
between the scales A 7 _i = l c x 2^ 1/2 and Aj = l c x 2 j+l/2 . The 
scales are the diameters of kernels used in calculating of the den- 
sity fields Dj-[ and Dj. 
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